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ABSTRACT 

We have surveyed four of the globular clusters in the Fornax dwarf galaxy (clusters 1, 2, 
3, and 5) for RR Lyrae stars, using archival F555W and F814W Hubble Space Telescope 
observations. We identify 197 new RR Lyrae stars in these four clusters, and 13 additional 
candidate horizontal branch variable stars. Although somewhat restricted by our short obser- 
vational baseline, we derive periods and light-curves for all of the stars in the sample, and 
calculate photometric parameters such as mean magnitudes and colours. This is the first time 
that RR Lyrae stars in the Fornax globular clusters have been quantitatively identified and 
measured. We find that the Fornax clusters have exceptionally large specific frequencies of 
RR Lyrae stars, in comparison with the galactic globular clusters. It is likely that Fornax 1 has 
the largest specific frequency measured in any globular cluster. In addition, the Fornax clus- 
ters are unusual in that their RR Lyrae populations possess mean characteristics intermediate 
between the two Oosterhoff groups defined by the galactic globular clusters. In this respect 
the RR Lyrae populations in the Fornax clusters most closely resemble the field populations 
in several dwarf galaxies. Fornax 5 has an unusually large fraction of RRc stars, and also 
possesses several strong RRe (second overtone pulsator) candidates. 

With a large sample of horizontal branch variable stars available to us, we revise previ- 
ous measurements of the horizontal branch morphology in each cluster. The Fornax clusters 
most closely resemble the "young" galactic halo population defined by Zinn in that their hor- 
izontal branch morphologies are systematically redder than many galactic clusters of similar 
metallicity. We also confirm the existence of the second parameter effect among the For- 
nax clusters, most markedly between clusters 1 and 3. The edges of the instability strip are 
well defined in several of the Fornax clusters, and we are able to make measurements of the 
intrinsic V — I colours of these edges. Finally, we determine foreground reddening and dis- 
tance estimates for each cluster. We find a mean distance modulus to the Fornax dwarf of 
(m — M)o = 20.66 ± 0.03 (random) ±0.15 (systematic). Our measurements are consistent 
with a line of sight depth of ^ 8 — 10 kpc for this galaxy, which is in accordance with its di- 
mensions as measured in the plane of the sky. This approximately spherical shape for Fornax 
is incompatible with tidal model explanations for the observed high internal stellar velocity 
dispersions in many dwarf spheroidal galaxies. Dark matter dominance is suggested. 

Key words: stars: variables: other - stars: horizontal branch - galaxies: star clusters - glob- 
ular clusters: general - galaxies: individual: Fornax dwarf spheroidal 



1 INTRODUCTION 

The Fornax dwarf galaxy is one of the most massive of the dwarf 
spheroidal (dSph) galaxies associated with the Milky Way (second 
only to the disrupted Sagittarius dSph), and the only undisturbed 
local galaxy to possess globular clusters (although Kleyna et al. 
(2003) have recently suggested the presence of a disrupted cluster 
in the Ursa Minor dSph). In this respect Fornax is somewhat un- 
usual, because its population of five means that it has the highest 
specific frequency of globular clusters for any known galaxy. 
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As summarized by Strader et al. (2003), although small in 
number, this globular cluster system is remarkably complex. Buo- 
nanno et al. (1998b; 1999) have shown that clusters 1, 2, 3, and 
5 appear coeval with both each other and the metal-poor galactic 
globular clusters, and that cluster 4 may be up to three gigayears 
younger. Strader et al. (2003) suggest that Fornax 5 may also be 
slightly younger. It has also been found (Smith et al. 1996; Buo- 
nanno et al. 1998b) that the Fomax globular clusters show evidence 
of the second parameter effect, whereby another parameter in addi- 
tion to metallicity appears to define the morphology of the horizon- 
tal branch in a cluster. It is often assumed that age is the second pa- 
rameter; however, if the age measurements listed above are correct, 
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then it is possible that age is not the sole second parameter in the 
Fornax clusters. We (Mackey & Gilmore 2003b) have studied the 
surface brightness profiles of the Fornax globular clusters, and find 
a large amount of variation in the cluster structures - ranging from 
an extremely extended cluster (Fornax 1), to a strong post core- 
collapse candidate (Fornax 5). Evidently the Fornax cluster system 
is highly complicated; however both this complexity, and the sys- 
tem's isolated nature provide unique opportunities for addressing 
the outstanding difficulties in our understanding of globular cluster 
formation and evolution. The Fornax globular clusters are therefore 
worthy of detailed study. 

Much can be learned about a stellar population from the vari- 
able stars within that population, and given the relatively close 
proximity of the Fornax dwarf, it is perhaps surprising that a de- 
tailed investigation of stellar variability in its clusters has not been 
undertaken. In fact, we can find no quantitative identification and 
measurements of any variable stars in these objects in the literature 
- only several passing mentions of possible RR Lyrae detections 
exist (Buonanno et al. 1985; Demers et al. 1990; Buonanno et al. 
1996; Smith et al. 1996; Smith et al. 1997; Buonanno et al. 1998b), 
along with one short article describing a survey in progress, and 
including cluster 3 (Maio et al. 2003). 

In the process of obtaining photometry for colour-magnitude 
diagrams from archival Hubble Space Telescope observations, we 
noticed the possibility of surveying four of the Fornax globular 
clusters for RR Lyrae stars. We present here the results of this sur- 
vey, which has been successful beyond our original expectations. 
We have discovered 197 RR Lyrae stars in the four clusters (For- 
nax 1, 2, 3, and 5), as well as 13 candidate horizontal branch vari- 
ables. We describe the data and our survey strategy in Section 2, 
and present light curves and mean properties for each star in Sec- 
tion 3. Finally, we have been able to use the RR Lyrae discoveries to 
measure detailed information about each of the four clusters under 
consideration. This discussion is presented in Section 4. 



2 OBSERVATIONS AND REDUCTIONS 
2.1 Data 

Wide Field Planetary Camera 2 (WFPC2) observations of four of 
the globular clusters in the Fornax dwarf galaxy (clusters 1, 2, 3, 
and 5, according to the notation of Hodge (1961)) are available 
in the HST archive (Program ID 5917). These observations were 
made between 1996 June 4 and 1996 June 6, through the F555W 
and F814W filters. The observational details are listed in Table 1. 

Each cluster was imaged 14 times in F555W and 16 times in 
F814W, with individual exposure durations ranging from 160 — 
600s, and 120 — 900s, respectively. Long total integration times 
are necessary because of the faintness of the Fornax globular clus- 
ters, which have distance moduli (m — M)o ~ 20.68 (Buonanno et 
al. 1999). The image sampling, presumably primarily intended to 
facilitate the removal of cosmic rays over the long integrations, ren- 
ders the data suitable for identifying sources with short period vari- 
ability, such as RR Lyrae stars. The exposure durations complement 
such identification, with horizontal branch (HB) stars being among 
the brightest stars in each image, and neither saturated nor overly 
faint in any image. Furthermore, observations through the two fil- 
ters are interleaved, so that for any given F555W observation, it is 
possible to find an F814W image with a similar exposure duration 
and an observation date matching to ±0.02 days. Such pairing of 
images allows colour information over a star's variability cycle to 
be measured. 



Although otherwise well suited for an RR Lyrae star search, 
the data have very short baselines of observation. This means that 
while variability can be easily detected, period determination is not 
straightforward - especially for RRab type variables, which typi- 
cally have periods in the range 0.4 — 1.1 days. Unless such stars 
are fortuitously measured at a critical region of their light curve, 
the under-sampling renders accurate period determination virtually 
impossible. The shorter cycle RRc type variables, which have peri- 
ods of 0.2 — 0.5 days, are more suited to the observation baselines. 
We will return to the question of period determination in Section 
2.4. 

Finally, we note that there also exist archival WFPC2 obser- 
vations of Fornax cluster 4 (Program ID 5637). However, the data 
consist of only three images per filter, and given the significant 
crowding and field star contamination for this cluster, were not suit- 
able for the present RR Lyrae work. We do not consider Fornax 4 
further. 

2.2 Photometry 

All WFPC2 images we retrieved from the HST archive under- 
went preliminary reduction according to a standard pipeline, us- 
ing the latest available calibrations. For photometric measurements 
from these calibrated images we used Dolphin's HSTphot (Dol- 
phin 2000a). This is a photometric package specifically designed 
for use on WFPC2 data. In particular, it uses a library of point 
spread functions (PSFs) tailored to account for the undersampled 
WFPC2 PSFs, allowing very accurate stellar centroiding and pho- 
tometry. 

Before running HSTphot, we completed several preprocessing 
steps. For each cluster, a reference image was identified (these are 
listed in Table 1) and all other images aligned with this image. This 
was achieved using the IRAF task IMALIGN, treating each image 
offset as a simple x- and y-shift. This is perfectly adequate for the 
present data - each image for a given cluster is offset from the 
reference image due to a simple dithering pattern of several pixels. 
There are no significant higher order distortions to account for. 

Next, we used the utility software accompanying HSTphot to 
mask image defects (such as bad pixels and columns), calculate a 
background sky image (for use in the photometry calculations), and 
attempt to remove cosmic rays and unmasked hot-pixels. We then 
made photometric measurements on each image using HSTphot in 
PSF fitting mode (as opposed to aperture photometry mode). We 
set a minimum threshold for object detection of 3u above the back- 
ground. We also enabled two additional features of HSTphot. First, 
we elected to calculate a local adjustment to the background image 
before each photometry measurement, using the pixels just beyond 
the photometry radius (Dolphin 2000a). This helped account for 
regions of rapidly varying background, such as near the centres of 
the very crowded clusters 3 and 5. Second, we chose to run artificial 
star tests in conjunction with the photometry measurements. This 
option generates a large number of "fake" stars on a CMD, places 
them one at a time on the image, and solves each just as if it were a 
real star. This allows quantitative estimates of detection complete- 
ness as a function of a star's magnitude, colour, and position in a 
cluster. 

For each photometric measurement of a star (ideally, 
14xF555W measurements and 16xF814W measurements) HST- 
phot output nine parameters, including flight magnitudes, standard 
magnitudes, and quantities characterizing the goodness-of-fit of the 
PSF solution (e.g., classification type, \, S/N, and sharpness). Mag- 
nitudes from HSTphot are calibrated according to the recipe of 
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Table 1. WFPC2 observations of four of the globular clusters in the Fornax dwarf galaxy (Program ID 5917). 



Cluster 


Data-group 


Reference 


Baseline 




Exposure 


Data-group 


N irn 


Exposure 


Name 


F555W 


Image 


(days) 




Durations (s) 


F814W 




Durations (s) 


Fornax 1 


u30m010eb 


u30m010et 


0.3595 


14 


3 x 600s, 4 x 500s, 


u30m010ib 


16 


2 x 900s, 6 x 700s, 












3x400s,4xl60s 






2 x 500s, 6 x 120s 


Fornax 2 


u30m020eb 


u30m020et 


0.3602 


14 


3 x 600s, 4 x 500s, 


u30m020ib 


16 


2 x 900s, 6 x 700s, 












3x400s,4xl60s 






2 X 500s, 6 X 120s 


Fornax 3 


u30m030eb 


u30m()30et 


0.3602 


14 


3 x 600s, 3 x 500s, 3 x 400s, 


u30m030ib 


16 


2 x 900s, 6 x 700s, 












Ix378s,4xl60s 






2 x 500s, 6 x 120s 


Fornax 5 


u30m040eb 


u30m040et 


0.3553 


14 


3 x 600s, 4 x 500s, 


u30m040ib 


16 


2 x 900s, 6 x 700s, 












3x400s,4xl60s 






2 x 500s, 6 x 120s 



Holtzman et al. (1995), and using the latest updates of the Dol- 
phin (2000b) CTE and zero-point calibrations. Each magnitude is 
corrected to a 0'.'5 aperture. HSTphot also provided mean posi- 
tional information for each object, in the form of a chip number 
(where chip refers to the PC, and chips 1 - 3 to WFC2-WFC4 
respectively) and pixel coordinates relative to the frame of the ref- 
erence image. These coordinates, in conjunction with the IRAF task 
METRIC and the cluster centres measured in Mackey & Gilmore 
(2003b), allowed us to calculate the radial distance of a given ob- 
ject from the centre of its cluster. 

We used the goodness-of-fit parameters to select only objects 
with high quality photometry. We retained only measurements for 
which an object was classified as stellar (HSTphot types 1, 2, and 
3), and for which % < 2.5, S/N > 3.0, -0.3 < sharpness < 0.3, 
and errors in the flight magnitude of < 0.1. After the application 
of this filter, we kept only stars with five or more F555W-F814W 
pairs of successful measurements (where a pair is defined as two 
observations within ±0.02 days of each other, as described in Sec- 
tion 2.1). We also passed the artificial star measurements through 
the quality filter. This allowed us to assign a detection complete- 
ness a c to each real star, by finding the fraction of successful arti- 
ficial star measurements in a brightness-colour-position bin about 
the real star. Bin widths were 0.25 mag in brightness, 0.2 mag in 
colour, and 2" in radial distance from the cluster centre. As an ex- 
ample, we consider a typical Fornax globular cluster RR Lyrae star, 
of magnitude V = 21.3, colour V — I = 0.5, and radial distance 
r — 7" . The fraction of successful artificial star measurements with 
21.175 < V < 21.425, 0.4 < V - I < 0.6, and 6" < r < 8" 
defines the detection completeness for this star. We expect the com- 
pleteness values so derived to be accurate to a few per cent, except 
for very low fractions. Values of a c less than ~ 0.25 should be 
regarded with caution. 



V 



2.3 Identification of Horizontal Branch variable stars 

The completion of the photometry procedure resulted in a list of 
stars in each cluster, each with between 5 and 14 pairs of F555W 
and F814W measurements. In general, most stars possessed ~ 14 
measurement pairs. Using colour-magnitude diagrams (CMDs) we 
selected stars in the HB region of each cluster, and determined the 
Welch & Stetson (1993) variability index Iws for each. This in- 
dex is calculated as follows. For N a epochs of observations, each 
resulting in measurements Vj ± av,i and U ± oi ti , the weighted 
mean magnitudes 



E«o 1 



r 



may be computed. The variability index is then defined as: 



where the normalized magnitude residuals are 



5Vi = 



Vi-V 

0~V : i 



5L = 



Ii-I 

07, z 



(1) 



(2) 



(3) 



The variability index therefore allows one to preferentially search 
for variations in photometric measurements which are correlated 
in the two colours being measured. For example, an RR Lyrae star 
will brighten and fade in both V and /, in synchronization over its 
pulsation cycle. Variability such as this results in a large value for 
Iws- In the case of random errors - for example where one or more 
stellar images has been impacted by a cosmic ray, causing spurious 
variability in one colour band, or in the case where crowding intro- 
duces large random errors into each measurement - then SVi and 
Sli are uncorrelated, and the expectation value of Iws is zero. 

We therefore selected stars for which Iws was greater than 
some positive constant as candidate variables. A small amount of 
experimentation showed Iws > 3.0 to be a suitable value. The 
resultant candidate variable stars are shown in Fig. 1. We identi- 
fied 283 candidates in total (24, 56, 137, and 66, in clusters 1, 2, 
3, and 5, respectively). In each of the HB regions on the CMDs, 
the RR Lyrae strip is evident, and in most cases clearly defined on 
both edges. For each cluster, there also in general exist some can- 
didate variable stars in the blue HB (BHB) region (except Fornax 
1, for which no BHB extension is observed), and some along the 
red giant branch (RGB). Because we did not want to assume any 
information, a priori, about the exact limits of the instability strip 
on the HB, we tested all flagged variable stars equally as RR Lyrae 
candidates. 



2.4 Light-curve fitting 

The best way to identify a candidate star as an RR Lyrae, and in ad- 
dition obtain quantitative information about its variability param- 
eters, is to attempt to fit a light curve to its photometry. For very 
under-sampled data, such as the present measurements, this can be 
tricky. Eschewing elegance, the most reliable method is the simple 
application of brute force. We used a routine written by Andrew 
Layden (described in Layden & Sarajedini (2000); Layden et al. 
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Figure 1. Colour-magnitude diagrams showing the horizontal branches of 
each of the four clusters. Stars with Iyys > 3.0 are marked as bold points 
- these are the candidate variables. The RR Lyrae region on each horizontal 
branch is well defined. We identified a total of 283 candidate variables - 24, 
56, 137, and 66, in the four clusters respectively. 



(1999); Layden (1998); and the references therein) for our period 
measurements. 

This program works by taking a set of 10 variable star tem- 
plates and attempting to fit them to the F-band photometric mea- 
surements for a number of different periods. The periods are se- 
lected by an incremental increase AP over a range Pi to P2. 
For each candidate period, the measurements are folded, and a 3- 
parameter fit (magnitude zero-point, phase zero-point, and light- 
curve amplitude) made to each template. Six of the ten templates 
represent RRab type variable stars, while two represent RRc type 
pulsators, and two represent variable binary stars (a W Ursae Ma- 
joris contact binary, and an Algol eclipsing binary). The RRab 
curves (templates 1 — 6), and one of the RRc curves (template 7) 
are derived from high quality measured light curves, as described in 
Layden & Sarajedini (2000), and Layden (1998). The second RRc 
template (number 8) is a simple cosine curve. From each template 
fit, x 2 1S calculated. Period-template combinations with small \ 2 
are likely to be the best representations of the measured data. 

Such a procedure works well when the sampling is spread over 
a much longer baseline than the period of the variable star under 
consideration. In the present case however, our ~ 0.36 day baseline 
is in general less than one RR Lyrae pulsation cycle. For the short 
period (RRc) variables, which typically have P < 0.45 days, this 
does not pose too much of a problem. It is possible to derive a 
reliable period estimate for these stars, although naturally, accuracy 
would be greatly increased by having a baseline of many periods. 

The situation is not so good for the RRab variables, which 
might have periods typically about twice that of our baseline. It 
is possible to determine a reliable period estimate for such stars, 
but only under certain circumstances. Specifically, it is necessary 
to observe the full amplitude of a star's variability (this allows the 
three parameters in a template fit to be well determined) and some 



portion of its decline from maximum (which allows discrimination 
between different potential periods). For a mean RRab period of 0.6 
days (in which case our baseline is 0.6P), examination of Figure 
1 in Layden (1998) shows that we would observe the appropriate 
portion of the light curve in approximately ~ 40 per cent of cases. 
For shorter period RRab stars this fraction is greater, and for longer 
period stars it is smaller. Another ~ 10 — 20 per cent of the time, 
we might observe most of the portion of the cycle which we re- 
quire, in which case the derived period will be a good estimate. In 
the remainder of cases, without assuming some prior knowledge of 
either the star's period or its amplitude, we cannot make an accurate 
period determination. In these incidences, we must make do with 
setting a lower limit to the period, which we achieve by assuming 
that the observed amplitude is the real amplitude of the star. The 
presence of these lower limits in our period measurements must 
be borne in mind if the data is put to any interpretive uses (see e.g., 
Section 4). Measurements at future epochs are required to pin down 
accurately the periods of such stars. 

Since we were searching for RR Lyrae stars, we chose a pe- 
riod range of 0.2 — 1.0 days, with AP = 0.01 per cent. For each 
candidate, we examined the few fits to the l^-band photometry with 
the smallest \ 2 ■ We used the J-band data as a consistency check, 
to verify that a fit which appeared suitable for the V-band pho- 
tometry also resulted in a good J-band light-curve. In general the 
/-band data was of significantly poorer quality than the V-band 
data, so fitting templates directly to these measurements did not add 
useful information to the procedure. In the vast majority of cases, 
the V-band solution with the very smallest x' 2 provided evidently 
the best fit. We conservatively estimate the typical random uncer- 
tainty in any of our fitted periods to be ~ 0.005 days. While we 
identified plenty of RR Lyrae stars, we did not find any stars with 
data clearly best fit by either of the variable binary star templates. 
However, given our restriction to a small portion of the CMD, and 
our very short observation baseline, this is not too surprising. Most 
of the candidate variables flagged on the BHB or RGB seemed to 
have either poor measurements or some form of irregular variabil- 
ity. Certainly none had values of Iws as large as those for typical 
RR Lyrae stars. 



3 RESULTS 

3.1 RR Lyrae stars 

Of the 283 candidate objects, we identified 197 as RR Lyrae stars 
(15, 43, 99, and 40, in clusters 1, 2, 3, and 5, respectively). For each 
of these stars, we obtained template fits and period measurements. 
The best fitting V-band light curves are presented in Fig. 2, and the 
positions, periods, and fit parameters (template, rms scatter) listed 
in Table 2. In this table, the coordinates (chip, X c , Y c ) are relative 
to the reference images listed in Table 1, and are in units of pixels. 
"Chip" refers to the camera on which the star was imaged - is 
the PC, and 1 - 3 are WFC2-WFC4. The radii (r) and complete- 
ness fractions (a c ) are calculated as described in Section 2.2. Nt 
refers to the label of the best fitting template - numbers 1 — 6 are 
RRab templates, and 7 — 8 are RRc templates. We consider periods 
marked with a star to be lower limit measurements. Note that on 
average, ~ 53 per cent of RRab stars have periods which are not 
lower limits - in accordance with our estimate in Section 2.4. In 
addition, we have marked several stars with question marks next to 
their type. In these cases, the light curves were best fit by a template 
of the RR Lyrae type listed (i.e., RRab or RRc), but could be fit by 
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Figure 2. V-band photometric measurements and best fitting light-curves for the 197 RR Lyrae stars. In addition the V-band photometric measurements 
for the 13 candidate variable stars are plotted. In each plot, the hollow circles represent points which have measurement errors more than 1.5 times larger 
than the mean photometric error for observations of the given star. These usually represent the shortest duration exposures in each set. Typical measurement 
errors are just within the plotted points. This photometry is available on-line at http : / /www . ast . cam. ac . uk/STELLARPOPS/Fornax_RRLyr/. 



a template of the other RR Lyrae type almost equally as well. Fur- 
ther measurements are required to confirm the listed classification 
for these stars. 

It is clear from the a c values that the only regions significantly 
affected by detection incompleteness are the very centres of the 
very crowded clusters 3 and 5. For cluster 3, completeness is sig- 
nificantly degraded within ~ 4", and there are only two detections 



within 2". Cluster 5 appears even more compact (cf. Mackey & 
Gilmore (2003b)), with completeness only degraded within ~ 2"; 
however, we note that there is only one detection within 1.5". 

From the best fitting light curves, we computed an intensity 
mean V magnitude, denoted (V), and the maximum, minimum, 
and amplitude of the variability (V+, and AV, respectively). 
We note that these quantities are in general not significantly dif- 
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Figure 3. Example finder charts for eight of the RR Lyrae stars - two per cluster. From left to right: (top) F1V02, F1V14, F2V21, F2V28; (bottom) F3V42, 
F4V49, F5V06, F5V18. Each chart is 15" on a side. Three of the stars (F1V02, F3V49, and F5V18) are on WFC chips; the remainder are on the PC. 



ferent from the same parameters derived directly from the F-band 
measurements. From the /-band measurements, we calculated the 
intensity mean (/}. iV is the number of pairs of observations for a 
given star. The parameter °E m is the epoch of maximum, in days 
since MJD 50 000. The reader should bear in mind that this param- 
eter is derived using a star's measured period. These quantities are 
all listed in Table 2. 

We also characterized the colour of each RR Lyrae star. In ac- 
cordance with Sandage (1990; 1993), we have calculated the mag- 
nitude mean of the N (V — I) measurements per RR Lyrae to 
best represent the colour of an "equivalent static star". This quan- 
tity, (V — I), is in preference to the two possible intensity mean 
colours (V) — (J) and (V — I). Finally, the magnitude mean colour 
of an RR Lyrae star at its minimum light (defined to be phases 
0.5 < (f> < 0.8) can be of use (see Section 4.4). For stars with 
measurements falling in this phase range, we also calculated this 
colour (denoted (V — I) m )- Both colours are listed in Table 2, to- 
gether with the calculated error a m in (V — I) m , and the number 
of observations at minimum light N m used in its calculation. 

Finally, using the coordinates listed in Table 2, we constructed 
a finder chart for each RR Lyrae star. For reasons of practicality, we 
have not included the full set of charts in this paper; however, some 
examples are displayed in Fig. 3. The remainder are available from 
the author on request. Alternatively, they may be downloaded from 
http://www. ast. cam. ac. uk/STELLARPOPS/Fornax_RRLyr/ along 
with the light-curves, and data from Table 2. 

3.2 Remaining candidates 

In addition to these 197 objects, a further 13 stars fell within or 
nearby to the HB instability strip defined by the RR Lyrae stars. 
However, suitable template fits could not be obtained for these^. 

t This was generally because the shape of the variation did not match any 
of the templates, or because a template which fit in the V-band produced a 



Plots of their photometry are presented at the end of Fig. 2. Most 
appear to have some form of cyclic variation, and given this and 
their locations on the HB, many are likely RR Lyrae stars of some 
sort. Their positions and completeness values are listed in Table 3, 
together with the photometric parameters (V), {/), V+, V_, AV, 
and N . Unlike for the RR Lyrae stars, the U-band parameters for 
the present stars were of course calculated from the measured data. 
Finder charts are also available for these candidate RR Lyrae stars 
on-line, or by request. 

3.3 Comparisons with previous work 

Although we searched carefully, we were unable to find any pub- 
lished quantitative measurements of RR Lyrae stars in Fornax glob- 
ular clusters. There has however been some work on the RR Lyrae 
population of the dwarf galaxy itself (e.g., Bersier & Wood (2002)). 
Nonetheless, several authors have noted the inevitable presence of 
RR Lyrae stars in the Fornax clusters, given their well populated 
horizontal branches, and on several occasions detections of vari- 
ability have been made. 

In one of the earliest CMD studies of the Fornax clusters, Buo- 
nanno et al. (1985) located several likely HB stars with variable 
photometry in their sample, and interpreted these as possible RR 
Lyrae stars. Buonanno et al. (1996) noted the presence of 17 hor- 
izontal branch stars with variable photometry in Fornax 1, and 40 
such stars in Fornax 3. They identified these as RR Lyrae stars; 
however there were not enough measurements to obtain significant 
average magnitudes. Smith et al. (1996) detected 21 candidate HB 
variables in Fornax 1 based on a comparison of their photometry 
with that of Demers, Kunkel & Grondin (1990), who also found 
evidence for variability in a couple of HB stars. Similarly, Smith, 

very poor light curve in the /-band. In many cases, the root cause is proba- 
bly poor measurements due to crowding or other contamination, combined 
with genuine variability 
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Table 2. Measured and calculated properties of the RR Lyrae stars. 
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0.482 


O O T A f\ f O 

237.4963 
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0.570 


RRab 


3 


0.030 


21.206 


20.687 
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0.421 


14 


21.040 


21.661 
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0.516 
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0.022 


2 


T^^l \ Tf\f\ 

F2V09 


1 


171.42 


249.32 


36.51 


1.00 
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0.612 


0.029 


5 


F2V11 





253.69 


501.65 


5.03 


1.00 
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0.387 


RRc 


7 
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0.014 


21.344 


20.655 
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0.98 
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0.044 


4 



Name 


Chip 


x c 

(pix) 


Y c 
(pix) 


r 

(") 




Period 
(days) 


Type 


N t 


rms 


(V) 


(I) 


F2V26 





614.26 


763.72 


18.65 


1.00 
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14 


20.709 


21.784 


1.075 


238.6336 


0.689 


0.016 


4 


F3V70 





207.17 


481.52 


8.38 


1.00 


0.664 


RRab 


5 


0.025 


21.308 


20.625 


0.539 


13 


20.867 


21.569 


0.702 


238.9016 


0.699 


0.031 


1 


F3V71 


2 


481.32 


510.67 


87.41 


1.00 


0.701 


RRab 


6 


0.031 


21.403 


20.852 


0.581 


14 


21.122 


21.612 


0.490 


238.4729 


0.642 


0.015 


8 


F3V72 





332.84 


570.22 


6.67 


0.89 


0.523* 


RRab 


6 


0.027 


21.055 


20.458 


0.590 


14 


20.817 


21.229 


0.412 


238.4555 


_ 


_ 





F3V73 





400.86 


340.76 


4.26 


1.00 


0.415* 


RRab 


6 


0.037 


21.337 


20.706 


0.580 


14 


20.929 


21.654 


0.725 


238.8752 


0.713 


0.036 


4 


F3V74 





220.40 


692.51 


13.93 


1.00 


0.412 


RRc 


7 


0.028 


21.266 


20.705 


0.564 


13 


21.046 


21.485 


0.439 


238.7453 


0.672 


0.022 


4 


F3V75 





382.96 


385.97 


2.14 


0.81 


0.391* 


RRab? 


6 


0.024 


21.142 


20.575 


0.542 


14 


20.962 


21.271 


0.309 


238.5093 


0.582 


0.035 


4 


F3V76 


2 


126.87 


746.74 


90.31 


1.00 


0.528 


RRab 


1 


0.049 


21.489 


20.819 


0.625 


14 


20.697 


21.980 


1.283 


238.3870 


0.694 


0.039 


2 


F3V77 





478.84 


443.13 


4.30 


1.00 


0.590 


RRab 


5 


0.025 


21.082 


20.476 


0.599 


14 


20.781 


21.252 


0.471 


238.4935 


0.632 


0.048 


5 


F3V78 





238.81 


426.41 


6.66 


0.91 


0.721 


RRab 


2 


0.020 


21.126 


20.549 


0.557 


14 


20.689 


21.434 


0.745 


239.0515 


0.646 


0.038 


4 


F3V79 


3 


364.46 


109.74 


49.93 


1.00 


0.655 


RRab 


3 


0.012 


21.485 


20.858 


0.572 


14 


21.117 


21.691 


0.574 


238.9895 


_ 


_ 





F3V80 


1 


276.06 


192.60 


32.34 


1.00 


0.621 


RRab 


1 


0.042 


21.265 


20.757 


0.624 


14 


20.602 


21.657 


1.055 


238.5658 


0.638 


0.015 


8 


F3V81 





53.44 


643.37 


17.76 


0.88 


0.657 


RRab 


4 


0.020 


21.488 


20.970 


0.504 


10 


21.006 


21.859 


0.853 


238.3446 


_ 


— 





F3V82 





443.58 


442.36 


2.70 


0.81 


0.509 


RRc 


7 


0.019 


21.046 


20.501 


0.500 


14 


20.935 


21.150 


0.215 


238.5657 


0.518 


0.024 


6 


F3V83 





295.07 


466.47 


4.37 


1.00 


0.444 


RRc 


7 


0.020 


21.197 


20.676 


0.545 


14 


20.975 


21.416 


0.441 


238.6740 


0.639 


0.019 


4 


F3V84 





419.13 


367.98 


3.34 


0.80 


0.312* 


RRab? 


1 


0.018 


21.345 


20.716 


0.596 


14 


21.128 


21.457 


0.329 


238.8683 


0.630 


0.034 


2 


F3V85 





437.33 


355.86 


4.24 


1.00 


0.509* 


RRab 


1 


0.040 


21.148 


20.691 


0.560 


12 


20.567 


21.483 


0.916 


238.6852 


0.632 


0.020 


3 


F3V86 


2 


539.43 


80.14 


69.02 


1.00 


0.648 


RRab 


1 


0.036 


21.504 


20.949 


0.632 


14 


20.942 


21.828 


0.886 


238.5342 


0.686 


0.023 


8 


F3V87 





358.82 


594.17 


7.42 


1.00 


0.397* 


RRab 


6 


0.017 


21.386 


20.753 


0.587 


14 


21.097 


21.601 


0.504 


238.4966 


0.666 


0.017 


4 


F3V88 


1 


647.16 


88.87 


47.12 


1.00 


0.524* 


RRab 


6 


0.029 


20.469 


19.968 


0.498 


14 


20.128 


20.728 


0.600 


238.8568 


0.595 


0.015 


4 


F3V89 





282.57 


458.47 


4.80 


1.00 


0.440* 


RRab 


6 


0.027 


21.379 


20.737 


0.616 


13 


21.107 


21.580 


0.473 


238.8521 


0.669 


0.021 


3 


F3V90 





434.05 


379.93 


3.29 


0.81 


0.664 


RRab 


1 


0.026 


21.214 


20.755 


0.597 


14 


20.593 


21.577 


0.984 


238.5428 


0.551 


0.022 


8 


F3V91 





432.48 


394.74 


2.78 


0.80 


0.453* 


RRab 


5 


0.033 


21.356 


20.727 


0.616 


14 


21.115 


21.488 


0.373 


238.8540 


0.637 


0.019 


4 


F3V92 





343.88 


428.18 


1.89 


0.04 


0.688 


RRab 


5 


0.027 


20.968 


20.408 


0.670 


14 


20.523 


21.231 


0.708 


238.5377 


0.705 


0.023 


7 


F3V93 


1 


645.48 


496.28 


74.83 


1.00 


0.536* 


RRab 


4 


0.026 


21.476 


20.851 


0.653 


14 


21.183 


21.688 


0.505 


238.7887 


0.670 


0.029 


4 


F3V94 





712.80 


327.71 


15.62 


1.00 


0.362* 


RRab? 


6 


0.011 


20.650 


20.061 


0.587 


14 


20.367 


20.861 


0.494 


238.8351 


0.677 


0.014 


3 


F3V95 





473.68 


389.27 


4.50 


0.88 


0.394* 


RRab 


6 


0.014 


21.090 


20.421 


0.652 


14 


20.907 


21.221 


0.314 


238.9062 


0.669 


0.016 


3 


F3V96 


1 


219.83 


284.31 


41.08 


1.00 


0.659 


RRab 


5 


0.029 


21.361 


20.816 


0.633 


14 


20.896 


21.638 


0.742 


238.5711 


0.629 


0.014 


5 


F3V97 





478.20 


456.30 


4.37 


0.82 


0.397* 


RRab 


6 


0.024 


21.361 


20.652 


0.661 


14 


21.110 


21.545 


0.435 


238.4847 


0.717 


0.028 


4 


F3V98 





552.32 


766.78 


16.89 


1.00 


0.451* 


RRab 


3 


0.057 


21.351 


20.669 


0.646 


14 


21.062 


21.508 


0.446 


238.4394 


0.692 


0.027 


4 


F3V99 





502.74 


389.31 


5.72 


0.82 


0.383* 


RRab 


6 


0.031 


21.368 


20.695 


0.649 


13 


21.083 


21.580 


0.497 


238.5170 


0.722 


0.019 


4 


F5V01 





420.01 


498.13 


2.69 


0.99 


0.405 


RRc 


8 


0.041 


21.007 


20.676 


0.345 


14 


20.771 


21.272 


0.501 


237.9404 


0.471 


0.025 


6 



Name 


Chip 


(pix) 


Y 

c 

Vpixj 


\ ) 




Period 


lype 


Nt 






IT\ 
\ J / 


F5V02 





91.38 


247.48 


16.23 


0.96 


0.311 


RRc 


7 


0.009 


21.493 


21.199 


F5V03 





445.78 


401.77 


3.43 


0.95 


0.278 


RRc 


7 


0.009 


21.326 


21.049 


F5V04 


3 


64.07 


311.29 


22.53 


0.91 


0.315 


RRc 


7 


0.019 


21.418 


21.041 


F5V05 





347.30 


687.40 


10.96 


1.00 


0.394 


RRc 


7 


0.015 


21.284 


20.920 


F5V06 





277.07 


471.31 


5.10 


0.98 


0.595 


RRab 


2 


0.011 


21.215 


20.631 


F5V07 





290.49 


593.95 


7.91 


1.00 


0.300 


RRc 


8 


0.017 


21.306 


20.904 


F5V08 





391.59 


375.00 


3.37 


0.99 


0.419 


RRc 


7 


0.023 


21.237 


20.848 


F5V09 





349.28 


424.61 


2.04 


0.99 


0.378 


RRc 


7 


0.014 


21.247 


20.836 


F5V10 





422.69 


317.87 


6.18 


1.00 


0.340 


RRc 


7 


0.020 


21.325 


20.915 


F5V11 





325.60 


391.00 


3.84 


0.99 


0.364 


RRc 


7 


0.018 


21.191 


20.786 


F5V12 


2 


323.81 


545.98 


81.04 


1.00 


0.396 


RRc 


7 


0.025 


21.143 


20.743 


F5V13 





396.96 


424.02 


1.22 


0.39 


0.368 


RRc 


7 


0.026 


21.273 


20.846 


F5V14 





380.64 


327.69 


5.53 


1.00 


0.372 


RRc 


7 


0.028 


21.317 


20.899 


F5V15 





465.75 


515.87 


4.70 


1.00 


0.362 


RRc 


8 


0.016 


21.291 


20.828 


F5V16 





405.52 


290.88 


7.24 


1.00 


0.520 


RRab 


3 


0.014 


21.332 


20.792 


F5V17 





554.79 


363.72 


8.56 


1.00 


0.353 


RRc 


7 


0.022 


21.310 


20.870 


F5V18 


1 


385.96 


81.77 


25.97 


1.00 


0.361 


RRc 


7 


0.015 


21.382 


20.960 


F5V19 


3 


378.12 


714.40 


71.86 


1.00 


0.404 


RRc 


7 


0.023 


21.47 


20.995 


F5V20 





694.33 


577.77 


15.10 


1.00 


0.422 


RRc 


7 


0.011 


21.356 


20.914 


F5V21 





417.89 


512.44 


3.21 


0.94 


0.351 


RRc 


7 


0.010 


21.171 


20.748 


F5V22 





468.89 


628.90 


8.97 


1.00 


0.402 


RRc 


7 


0.016 


21.367 


20.905 


F5V23 


1 


476.49 


152.30 


37.03 


1.00 


0.379 


RRc 


7 


0.033 


21.362 


20.91 


F5V24 





437.65 


424.47 


2.56 


0.99 


0.445 


RRc 


7 


0.015 


21.413 


20.934 


F5V25 





463.33 


340.50 


6.04 


1.00 


0.494 


RRab 


2 


0.020 


21.18 


20.685 


F5V26 





227.21 


589.26 


9.66 


0.98 


0.571 


RRab 


2 


0.016 


21.28 


20.712 


F5V27 





432.85 


616.59 


7.89 


1.00 


0.382 


RRc 


7 


0.014 


21.281 


20.837 


F5V28 





115.32 


737.13 


17.89 


0.95 


0.326 


RRc 


7 


0.021 


21.434 


21.073 


F5V29 





471.83 


291.22 


8.15 


1.00 


0.420 


RRc 


8 


0.019 


21.367 


20.874 


F5V30 


3 


211.47 


345.73 


37.21 


1.00 


0.433* 


RRab 


4 


0.019 


21.297 


20.728 


F5V31 





437.38 


321.31 


6.25 


1.00 


0.572 


RRab 


3 


0.012 


21.31 


20.706 


F5V32 





347.33 


243.78 


9.51 


1.00 


0.603 


RRab 


3 


0.032 


21.296 


20.794 


F5V33 





293.43 


457.39 


4.28 


1.00 


0.428 


RRc? 


7 


0.011 


21.234 


20.685 


F5V34 





326.89 


674.77 


10.60 


1.00 


0.720 


RRab 


3 


0.049 


21.433 


20.74 


F5V35 





351.30 


494.77 


2.64 


0.94 


0.415 


RRc? 


7 


0.014 


21.066 


20.547 


F5V36 


3 


501.81 


655.42 


77.55 


1.00 


0.540 


RRab 


3 


0.016 


21.488 


20.804 


F5V37 





378.58 


336.97 


5.11 


0.98 


0.532* 


RRab 


5 


0.015 


21.441 


20.923 


F5V38 





336.72 


619.27 


8.06 


0.98 


0.392* 


RRab 


5 


0.015 


21.524 


20.86 


F5V39 





616.54 


781.25 


18.24 


1.00 


0.513* 


RRab 


5 


0.017 


21.378 


20.836 


F5V40 





418.38 


479.85 


2.00 


0.95 


0.466* 


RRab 


6 


0.021 


21.373 


20.746 



(V-I) 


N 


v+ 


V- 






(V-I)m 


CTm 


0.343 


14 


21. 179 


o 1 om 
21.819 


0.640 


on nooo 
23 /.yoo2 


0.433 


0.032 


0.320 


14 


21.121 


21.528 


0.407 


OQO nCAQ 

238.0503 


0.410 


0.027 


n no 
O.J /o 


I4 


21.284 


O 1 C/1C 

21.343 


n oa 1 
0.201 


1 1 1 C A 

238.1334 


n a 1 1 
0.413 


n n 1 n 
0.01 / 


0.367 


13 


21.21 1 


21.351 


0.140 


OQO (\A1l\ 

238.0430 


0.384 


0.017 


0.536 


14 


20.713 


21.577 


0.864 


on "7i i 1 

237.7311 






0.377 


14 


21.205 


21.412 


0.207 


151 . 9470 


0.329 


0.020 


0.410 


14 


20.91 1 


21.578 


0.667 


no ft o i A 

238.0830 


0.532 


0.026 


0.400 


13 


21.016 


21 .477 


0.461 


oio nnAi 

2Jo.uyoi 


0.545 


0.025 


0.401 


14 


21.038 


21.619 


0.581 


238.137 / 


0.536 


0.026 


0.451 


14 


20.919 


21.468 


0.549 


o "5 1 ncn 

237.951 / 


0.548 


0.023 


f\ /1AC 

U.403 


1 A 

14 


on OIO 
20.835 


o 1 A cn 
21.439 


n ao 1 
0.021 


oio nr\Act 
238.0049 


n c c i 
0.331 


n no 1 
0.02 1 


0.439 


1 A 

14 


O 1 f\A A 

1 1 .U44 


o i cm 
21.301 


n A C*7 

0.43 / 


on n i nn 
23 / .9100 


n enn 
0.309 


n no h 
0.02 / 


n /iAn 

0.409 


1 A 

14 


o 1 n^? 
21.033 


O 1 COI 

21.383 


n coo 
0.328 


oio n/ino 
238.0498 


0.318 


n non 
0.020 


0.425 


14 


21.058 


21.552 


0.494 


on oTon 

237.8789 


0.568 


0.017 


U.430 


1 A 

14 


on nt A 
20. /10 


21 .708 


n nno 
0.992 


oio om 1 
2J0.2U1 1 


U.334 


n nn 
0.01 / 


0.397 


14 


21.055 


21.566 


0.5 1 1 


on onoT 

23 /.8923 


0.544 


0.034 


0.475 


14 


21.131 


21.635 


0.504 


oi ~7 nm i 

237.9931 


0.558 


0.019 


0.486 


14 


21.256 


21 .681 


0.425 


oio fiAin 
235.U03U 


0.570 


0.018 


0.416 


14 


21.068 


21.652 


0.584 


no ^Acn 

238.2059 


0.573 


0.019 


n a qo 
0.482 


1 A 

14 


on cti 1 
20.931 


O 1 A 1 O 

21.412 


n aq 1 
0.481 


on nAn/i 
23 /.y094 


n ceo 
0.332 


n noc 
0.023 


0.482 


14 


21.129 


21.605 


0.476 


on n/1 a c 

237.9445 


0.577 


0.022 


n a cn 
0.439 


1 A 

14 


1 1 1 OA 

21.120 


o 1 enn 
21.399 


n a ni 
0.4/3 


oqo 1 
2io. 142 / 


n /iTA 
0.4/0 


n nn 
0.01 / 


n A ?o 
0.432 


1 A 

14 


O 1 0"70 

21.2 /o 


O 1 C/11 

21.343 


n oac 
0.203 


on oicn 
23 /.8330 


ft A CA 

0.430 


n ni /i 
0.034 


0.521 


14 


20.630 


21.584 


0.954 


On OAA/1 

23 7.5604 


0.517 


0.050 


n ci o 
0.338 


1 A 

14 


on £on 
20.080 


O l "70A 

21 . /20 


1 C\A A 

1 .040 


on "7om 
23 /. /202 






0.489 


14 


21.055 


21.506 


0.451 


oq o m n/i 

238.0194 


0.518 


0.021 


0.454 


13 


21.152 


21.722 


0.570 


o"> o n 1 i o 

238.0419 


0.471 


0.020 


0.501 


13 


O 1 O 1 o 

21.212 


21.534 


0.322 


on nim 

23 / .yiu/ 


0.602 


n nio 
0.038 


0.518 


14 


on "7i n 

20.730 


21.746 


1.016 


OQO OQA*7 

238.230/ 


0.644 


0.015 


n CC7 
0.3 3 / 


1 A 

14 


on QA 1 
20.841 


O 1 coo 
21.382 


n n A 1 
0. /41 


on An nn 

23 / .oyoo 


n aoo 
0.082 


n non 
0.029 


0.534 


14 


20.850 


21.553 


0.703 


238.1314 


0.573 


0.016 


0.538 


14 


21.125 


21.337 


0.212 


238.2763 


0.591 


0.021 


0.490 


13 


20.747 


21.864 


1.117 


238.2507 






0.536 


14 


20.951 


21.174 


0.223 


238.0014 


0.539 


0.023 


0.559 


14 


20.971 


21.792 


0.821 


237.6789 


0.674 


0.023 


0.551 


14 


21.009 


21.696 


0.687 


238.1578 


0.640 


0.015 


0.627 


14 


21.279 


21.659 


0.380 


237.8331 


0.656 


0.017 


0.542 


14 


21.168 


21.493 


0.325 


238.1786 


0.636 


0.031 


0.617 


14 


21.163 


21.527 


0.364 


238.1904 


0.627 


0.017 
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Table 3. Photometric measurements for the candidate RR Lyrae stars. 



Name 


Chip 


X c 


Y c 


r 


Ctc 


(V) 


(-0 


No 


V+ 


V- 


AV 


F1C01 


2 


216.99 


73.27 


39.27 


1.00 


21.356 


20.796 


13 


21.017 


21.510 


0.493 


F2C01 





270.89 


750.76 


14.21 


0.98 


21.253 


20.729 


14 


20.995 


21.514 


0.519 


F2C02 





181.01 


437.51 


7.69 


0.97 


21.395 


20.877 


12 


21.286 


21.533 


0.247 


F2C03 





494.90 


97.50 


17.17 


1.00 


21.164 


20.577 


13 


20.924 


21.482 


0.558 


r jl,ui 


u 


OQO 1 ^ 
ZVZ. 1 j 


70 
J4J. 1 Z 


0.04 


1 .uu 


01 1/11 
Zl . 14 J 


on eno 


\A 
14 


01 f"K7 
Z 1 A) J 1 


01 0^^ 

zi .zoo 


u.zuv 


F3C02 





102.44 


473.52 


12.95 


1.00 


21.318 


20.825 


14 


21.069 


21.625 


0.556 


F3C03 





355.77 


453.05 


1.61 


0.08 


21.326 


20.851 


14 


20.832 


21.823 


0.991 


F3C04 





425.88 


793.12 


16.40 


0.88 


21.316 


20.714 
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Rich & Neill (1997) found evidence for 7 HB variables in Fornax 
3, and one variable above the Fornax 3 HB. Using the WFPC2 data 
from the present study, Buonanno et al. (1998b) identified 8, 39, 
66, and 36 candidate RR Lyrae stars in clusters 1, 2, 3, and 5, re- 
spectively, based on the fact that these stars had frame-to-frame 
variations greater than 3a. These numbers are quite similar to ours, 
as might be expected given that we are using the same observations. 
However, we detect more RR Lyraes in all clusters, demonstrating 
the sensitivity of our search strategy. 

More recently, Maio et al. (2003) have described wide-field 
observations of a region of the Fornax dwarf, including cluster 3. 
They detect 70 candidate variables in a 95" x 95" box centred 
on this cluster, most of which lie on the HB. This is again fewer 
variables than we detect in Fornax 3; however, the present WFPC2 
observations have superior resolution to the terrestrial wide-field 
observations. Maio et al. (2003) present preliminary uncalibrated 
light curves for two of their candidate cluster 3 RR Lyraes; however 
neither these nor any of their other variables are identified, and no 
quantification of the variability is presented. 

We therefore believe that all of the 197 RR Lyrae stars and 
13 candidates we describe here have been newly identified. This is 
the first time that calibrated light curves and variability data have 
been published for any RR Lyrae stars in the globular clusters in 
the Fornax dwarf galaxy. 

3.4 Effects of under-sampling 

We consider briefly here the effects that our short baselines of ob- 
servation are expected to have on the RR Lyrae parameters derived 
from the fitted light curves. As described in Sections 2.4 and 3.1, 
the primary effect is that ~ 50 per cent of the calculated RRab peri- 
ods are lower limits. This has knock-on effects for most of the other 
quantities computed for these stars, which must be accounted for if 
any use is to be made of the measurements (e.g., in determining 
distance moduli). 

For example, suppose an RRab star has a true amplitude 
(AV)t = {V-)t — (V+)t and a true period of pulsation P t . How- 
ever, because of our limited baseline of observation (< P t ), we 
have only measured AV < (AV) t for this star, and hence deter- 
mined a period P < P t (as described in Section 2.4). This will lead 
to a systematic error in our calculation of (V) from the fitted light 
curve. It is possible that we have sampled only the brightest portion 
of the light curve, so that our measurement of V+ is accurate, but 



in this case our assumed V- will be brighter than (V-)t and our 
calculated (V) will also be brighter than the true intensity-mean V 
magnitude (V) t . Similarly, it is possible that we have accurately 
measured V-, but assumed too faint a value for V+. In this case our 
(V) will be fainter than (V)t- For symmetric light curves, such sys- 
tematic errors should average to zero over a large ensemble of stars; 
however RRab stars have significantly asymmetric pulsations, in 
that they are brighter than (V)t for typically only ~ 30 — 40 per 
cent of a cycle. Hence, while determinations of (V) for each star 
individually may be either too bright or too faint, we expect on av- 
erage that our measurements will provide values of (V) somewhat 
fainter than (V) t (i.e., we measure V- accurately more often than 
we measure V+ accurately). The same reasoning applies for calcu- 
lations involving (J), which are made directly from the photometry. 

For stars with lower limit periods, our measurement of the 
magnitude-mean colour over the cycle will also suffer a system- 
atic error. Because we will have preferentially sampled the RRab 
light cur ves around minimum light, on average a measurement of 
(V — I) will be redder than its true value, because minimum light 
is the reddest part of the pulsation cycle. We note that individual 
measurements may, of course, also be too blue - this is just more 
unlikely than a red error. 

It is more complicated to estimate the effect of a lower-limit 
period on a measurement of (V — I) m . Only points with phases 
in the range 0.5 < 4> < 0.8 are used in the calculation of this 
parameter. Since we have P < Pt, the true phase of any given point 
will be smaller than the phase calculated using P. Hence, it is likely 
that we have used points which have true phases <j> t < 0.5 in our 
calculations, and that on average our measured (V — I) m values 
will be too blue. We note that this effect will be small, since the 
rate of change of colour around the minimum of an RRab pulsation 
cycle is slow. 

There are two ways to account for these systematic errors if 
any use of the derived data is made - for example, in determining 
the distance modulus to Fornax. The first is simply to leave the 
stars flagged as having lower limit periods out of the calculations. 
Alternatively, it is possible just to use those stars which have large 
amplitudes (say AV > 0.9), since these are the most likely to 
have an accurate determination for P. These two methods allow 
a consistency check to be made, and the magnitude of any of the 
systematic errors described above to be investigated. 
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4 DISCUSSION 

Because the RR Lyrae populations of the Fornax globular clusters 
have not been previously investigated in a quantitative manner, it is 
useful to consider their characteristics. Much can be learned from 
the measurements we have made, even bearing in mind the uncer- 
tainties in the period fitting. As discussed in Section 3.4, it is pos- 
sible to investigate and account for the systematic errors our proce- 
dure might introduce. Throughout the following Sections, Table 4 
keeps a record of the calculations we discuss. 

4.1 Specific frequencies 

Each of the four Fornax clusters we have examined has a significant 
population of RR Lyrae stars. To investigate this quantitatively we 
calculate the RR Lyrae specific frequency, S RR , which is the num- 
ber of RR Lyrae stars in a cluster normalized to a cluster absolute 
magnitude of Mv = —7.5: 

S RR = N RR W 0A< - 7 - 5+M ^ (4) 

Buonanno et al. (1998b) have noted that S RR for Fornax clusters 
1, 2, 3, and 5 appears quite high relative to the Milky Way globular 
clusters. 

We calculate S RR using the integrated luminosities for the 
four Fornax clusters derived by us from accurate surface bright- 
ness profiles (Mackey & Gilmore 2003b). For the four clusters 
we obtained integrated luminosities log(L<x> / '£©) = 4.07, 4.76, 
5.06, and 4.76 respectively. These correspond to F555W luminosi- 
ties. Assuming the F555W solar magnitude to be Vgfg = +4.85 
(Mackey & Gilmore 2003a), we obtain absolute magnitudes of 
—5.33, —7.05, —7.80, and —7.05 for the four clusters. These are in 
good agreement with the values derived by Buonanno et al. (1985), 
and the values listed by Webbink (1985). 

The specific frequencies for the four clusters are then S RR = 
110.7, 65.1, 75.1, and 60.5, respectively. We note that these are 
lower limits, because we have not included information about de- 
tection completeness, or counted the candidate variable stars. We 
can account for completeness effects using the a c values in Ta- 
ble 2. To avoid unwarranted weighting by the few stars with very 
low completeness fractions, we set any a c values of less than 0.25 
equal to 0.25. If we then assume that all the candidate variable 
stars are actually RR Lyrae stars, we obtain specific frequencies of 
S RR = 118.2, 70.5, 90.2, and 68.4. These are still lower limits, 
because in no cases did the WFPC2 observations cover an entire 
cluster. Hence, the samples of RR Lyrae stars are certainly incom- 
plete, but in calculating S RR , we have normalized the N RR values 
to full cluster luminosities. Accounting for this effect is not trivial, 
because we do not know the spatial distribution of RR Lyrae stars 
in any cluster, and because the WFPC2 geometry is complicated. 
Given that we have imaged the centre of each cluster, we expect 
the N RR to be 80 — 90 per cent complete. Hence, the S RR values 
may be 10 — 20 per cent greater than quoted above. 

These specific frequencies are extremely high compared to 
those for the galactic globular clusters. The catalogue of Harris 
(1996) lists only four of the ~ 150 galactic clusters as having 
S RR > 60, and only two of these have S RR > 100. The largest 
value, that for Palomar 13, is S RR = 127.5. Given our sample in- 
completeness, it is very likely that Fornax 1 has S RR significantly 
greater than this. It is certainly intriguing that only a tiny fraction 
of galactic globular clusters have very high S RR , whereas all of the 
Fornax clusters we have investigated do - this is a fact worthy of 
further attention. 



4.2 Oosterhoff classification 

It is also instructive to investigate the relationship between the de- 
rived periods and amplitudes for our RR Lyrae populations. This 
can be done by plotting a period-amplitude diagram for each clus- 
ter. For galactic globular clusters, such diagrams fall naturally into 
two distinct groups - the Oosterhoff type-I and type-II clusters (Ool 
and OoII, respectively). Each group has characteristic properties, 
summarized in Smith (1995). Furthermore, there is in general a 
close relationship between period and amplitude for the RRab stars 
in a cluster - namely that the shortest period RRab stars have the 
largest amplitudes, while those with the smallest amplitudes have 
the longest periods. 

Period-amplitude diagrams for the four Fornax clusters are 
presented in Fig. 4. RRab stars with good periods are plotted as 
filled circles, while those with lower limit periods are plotted as 
crosses. RRc stars are plotted as triangles - those with uncertain 
classification in Table 2 are hollow points, while the remainder are 
filled. From examination of these diagrams, evidence of the lower 
limit period problem is immediately clear. The crosses form a dis- 
tinct population between the good RRab stars and the RRc stars 
in each diagram. Once good periods are determined for these stars, 
they will move into the regions occupied by the good RRab stars. In 
the best populated clusters, such as Fornax 3, evidence of the upper 
part of the RRab sequence can be seen, although the relationships 
certainly have significant scatter. Improved period measurements 
should tighten these sequences. 

From the period-amplitude diagrams it is not possible to de- 
termine conclusively whether the Fornax clusters resemble either 
the Ool or OoII clusters, or (like the old LMC clusters and several 
dwarf spheroidal galaxies) have properties which lie somewhere 
in between the two (e.g., Pritzl et al. (2002)). We can shed some 
light on this by calculating selected parameters for each of our RR 
Lyrae populations. Smith (1995) lists the Ool clusters as having 
(P ab ) = 0.55 d, (P c ) = 0.32 d, and N C /N RR = 0.17, and 
the OoII clusters as having (P ab ) = 0.64 d, (P c ) = 0.37 d, and 
N C /N RR = 0.44. Pritzl lists these parameters for several old LMC 
clusters. These have (P ab ) = 0.56 - 0.60 d, (P c ) = 0.33 - 0.38 
d, and N C /N RR = 0.3 — 0.5, placing them in general somewhere 
between the Ool and OoII groups. 

We calculated these parameters for the four Fornax clusters 
- the results are listed in Table 4. Our (P c ) values cover the range 
0.37—0.43 d, which are similar to those for OoII clusters. However, 
excepting Fornax 5, the fraction N C /N RR ranges from 0.3 — 0.38 
for the Fornax clusters, which is more similar to the intermediate 
LMC clusters than either of the Oosterhoff groups. Given the pres- 
ence of the many RRab stars with lower limit periods, it does not 
make sense for us to calculate (P a b) using all our RRab stars^ . It 
is far more useful to calculate this quantity using only those RRab 
stars with good periods. If we do this, we obtain values in the range 
0.57 — 0.61, again most similar to the intermediate LMC clusters. 
We must bear in mind however that these averages do not include 
many stars with long periods, since we have only measured lower 
limit periods for these stars. Hence, these mean values will likely 
increase once all stars have good periods. 

One final discriminator we can use is the transition period 
Ptr of RRab stars - that is, the shortest period that this class of 
stars have in a given cluster, without being RRc stars. Given the 
inverse period-amplitude relation for RRab stars, those with the 



t although we list this value in Table 4 for comparative purposes. 



© 2003 RAS, MNRAS 000, 000-000 



20 A. D. Mackey & G. F. Gilmore 



> 
< 



3 



0.5 




0.4 0.5 
Period (days) 




0.4 0.5 0.6 0.7 

Period (days) 



1 - 



3 



E 



0.5 



1 1 1 

Fornax 


i i i i i i i i i 1 1 1 1 1 1 1 M 

3 

• 




• 

• x • • 

x 


A * * 


'• :• . 

x* x X • •• - 

A x X • 

A 

xfX " X 5 X A * • 


A A 


A X • 

*»x* * • 

A * X 


A i 

i , 


A \ X* • 

a * a 

A 

1 1 1 1 1 1 1 , 1 1 1 1 1 1 1 1 1 1 1 1 1 



0.3 0.4 0.5 0.6 0.7 

Period (days) 



P 



Q, 

E 



0.5 



1111 



Fornax 5 



^4 ' 



0.3 0.4 0.5 0.6 0.7 

Period (days) 



Figure 4. Period-amplitude diagrams for each of the four clusters. RRab stars with good periods are plotted as solid circles; those with lower limit periods 
are crosses. RRc stars with uncertain classification are hollow triangles, while the remainder are solid triangles. 



shortest periods have the largest amplitudes. Calculating (P a b) for 
our RRab stars with AP > 0.9 therefore provides a good esti- 
mate of P tr . Smith (1995) lists P tr = 0.43 for Ool clusters and 
Ptr = 0.55 for OoII clusters. Our values range from 0.54 — 0.58, 
which match the OoII clusters well. 

Hence we can rule out the possibility that the Fornax globu- 
lar clusters 1, 2, 3, and 5 resemble the galactic Ool type clusters. 
However, with quite a dispersion in their properties, it is not clear 
whether they do resemble either the OoII type clusters, or the in- 
termediate LMC type clusters. Both N c /Nrr and (P a b) suggest 
an intermediate-type classification, while (P c ) and P tr suggest a 
similarity to the OoII clusters. One explanation is that the Fornax 
clusters are, like the LMC clusters, of intermediate classification, 
but with somewhat different characteristics. This interpretation is 
consistent with the observation that the properties of RR Lyrae stars 
measured in the fields of dwarf galaxies classify these objects also 



as intermediate, but with a considerable range in properties (see 
e.g., Pritzl et al. (2002)). Bersier & Wood (2002) have measured 
the properties of the field RR Lyrae population in Fornax, and find 
(Pab) = 0.585 and a mean RR Lyrae [Fe/H] ~ — 1.8§. This metal- 
licity is very similar to those for Fornax clusters 2 and 5 (see Sec- 
tion 4.6 and Table 4), which also have (P ab ) ~ 0.58. Clusters 1 
and 3 are somewhat more metal poor, with [Fe/H] ~ — 2, and have 
slightly longer mean RRab periods - (P a b) ~ 0.61. From the list in 
Pritzl, these are most similar to the Carina and Draco dwarf galax- 
ies, which also have [Fe/H] ~ —2, and (P a b) ~ 0.62. In fact, 

§ In fact, they measure [Fe/H] ~ —1.6 on the Butler-Blanco metallicity 
scale, which is ~ 0.2 dex more metal rich than the Zinn & West (1984) 
scale. We quote the metallicity on the latter scale for consistency with the 
measurements in Table 4. 
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the Carina RR Lyrae stars have recently been studied in detail by 
Dall'Ora et al. (2003), who find that (P a b) for these stars matches 
well the OoII clusters, but that N c /Nrr matches the Ool clusters 
- a similar situation to that for the Fornax clusters. These combina- 
tions are the strongest indicators that the Fornax globular clusters 
also deserve an intermediate Oosterhoff classification. 



4.3 The strange case of Fornax 5 

It is worth commenting briefly on cluster 5, which has somewhat 
exceptional properties. Specifically, it has a very large relative pop- 
ulation of RRc stars, so that N c /Nrr = 0.68. This is significantly 
larger than any of the RR Lyrae populations listed in the compila- 
tion of Pritzl et al. (2002), in which the largest belongs to the LMC 
cluster NGC 2257, with N C /N RR = 0.50. 

The period-amplitude diagram for Fornax 5 is also unusual, 
in that in addition to the well defined RRab and RRc groups, there 
appears a third group of four stars in the very lower left of the di- 
agram. This group strongly resembles the putative RRe type stars 
(second overtone pulsators) identified by Clement & Rowe (2000) 
in uj Centauri. Although there has been some debate in the liter- 
ature over the existence (or not) of such objects, there have been 
several instances of observational evidence (see the references in 
Clement & Rowe). The ~ 20 RRe candidates identified by Clement 
& Rowe have a mean period of 0.304 days, and a range 0.28 — 0.35 
days. The four candidates in Fornax 5 (F0503, F0504, F0505, and 
F0507) have a range 0.28 - 0.39 days, and a mean of 0.32 days. 
The mean amplitude of the uj Cen objects is slightly larger than 0.2 
mag, which also matches the mean amplitude of the Fornax 5 can- 
didates: (AV) ~ 0.25 mag. Hence the two groups of stars occupy 
the same region on the period-amplitude diagram, and we conclude 
that the four Fornax 5 objects are good RRe candidates. We note 
that these four stars are among the five bluest RR Lyrae stars in 
this cluster, adding weight to the argument. There also exist several 
similar stars in clusters 2 and 3; however these do not seem as well 
separated on the period-amplitude diagrams as the cluster 5 stars, 
and we cannot claim them as RRe candidates with certainty. 

Returning to the question of the large number of RRc stars 
in Fornax 5, even subtracting the four potential RRe stars from 
the calculation leaves N c /Nrr = 0.58 - still an unusually large 
value. It is interesting to note that Fornax 5 is unique in additional 
ways - it is in the very outer regions of the Fornax dwarf, and it 
is the only post core-collapse candidate in the system (Mackey & 
Gilmore (2003b)). It seems difficult to imagine some link between 
the structure of the cluster and its RR Lyrae star population how- 
ever, especially given that neither uj Cen nor NGC 2257 are excep- 
tionally compact objects. 

4.4 Reddening towards the Fornax globular clusters 

It is a well known property of ab-type RR Lyrae stars that they have 
extremely similar intrinsic colours at minimum light, independent 
of metallicity: [(V - I) m ]o = 0.58 ± 0.03 (see e.g., Mateo et 
al. (1995) and references therein). This can be used to determine 
the line-of-sight reddening to a group of RRab stars. The validity 
of this technique has been demonstrated for RR Lyrae stars in the 
Sagittarius dwarf galaxy by both Mateo et al. (1995) and Layden & 
Sarajedini (2000). 

To apply this technique to the four Fornax clusters individu- 
ally, we calculated the weighted mean {(V — I) m ) for each sample 
of RRab stars, using the (V — I) m and a m measurements from 



Table 2. We obtained ((V - I) m ) = 0.68 ± 0.01 for Fornax 1, 
0.65±0.01 for Fornax 2, 0.63±0.01 for Fornax 3, and0.62±0.01 
for Fornax 5. To investigate the possibility of systematic errors due 
to the lower limit periods, as described in Section 3.4, we recalcu- 
lated these values - first using only RRab stars not flagged in Table 
2 as having lower limit periods, and second using only RRab stars 
with AV > 0.9. The results are listed in Table 4. As expected, 
for the most part the values calculated using all the RRab stars are 
slightly bluer than those which exclude the stars with lower limit 
periods. However, the offsets are not constant. We note that the re- 
calculations for Fornax 1 and Fornax 5 each involved less than 5 
stars (since several objects do not have (V — I) m measurements), 
so these two values are strongly subject to stochastic effects. Over 
the combined ensemble of all stars, the two recalculation values are 
8 x 10~ 4 mag redder and 3 x 10 ~ 3 mag redder than the original 
calculation, respectively. We therefore conclude that the systematic 
effect is small (as predicted in Section 3.4), and we feel confident 
in using the four original values of ((V — I) m ) quoted above. 

It is then a simple matter to show that for the four clusters 
E{V - I) = 0.10, 0.07, 0.05, and 0.04. Errors in each of these 
values are ±0.01 random and ±0.03 systematic (from the value 
of [(V — J) m ]o determined by Mateo et al. (1995)). Adopting the 
reddening laws calculated in Mackey & Gilmore (2003b), we used 
these colour excesses to calculate Ay and E(B — V) for each clus- 
ter. All of these quantities are listed in Table 4. The four E(V — I) 
values we obtained match previous measurements well. For exam- 
ple, Buonanno et al. (1998b) determined E(V - I) = 0.05, 0.09, 
0.05, and 0.08 (all ±0.06) for the four clusters, using their red giant 
branches. 

4.5 Horizontal Branch morphology 

Quantitative measurements of the numbers and colours of RR Lyrae 
stars in the Fornax clusters allows a new determination of their hor- 
izontal branch morphologies and types. Fig. 5 shows a CMD of the 
HB region in each cluster. The RR Lyrae stars and candidate vari- 
ables are plotted according to their identification - RRab stars are 
circles, RRc stars are triangles, and candidates are crosses. Stars 
with uncertain classification have hollow points (circles or trian- 
gles). For each cluster, the edges of the instability strip appear well 
defined. RRc stars fall toward the blue edge and RRab stars toward 
the red edge, although there appears to be some mixing of the two. 
This may be intrinsic, or it might be due to colour and classification 
uncertainties (although these effects should be small). 

We note that Fornax 3 has four RR Lyrae type objects which 
are ~ 0.5 mag brighter than the main HB. The identity of these 
objects is not clear. It might be that they are interlopers from the 
field of the Fornax dwarf itself - cluster 3 has the highest level of 
background stars of the four clusters considered here (Demers, Ir- 
win & Kunkel (1994)). Alternatively, they might be unusual cluster 
members, such as Population II cepheids or anomalous cepheids. 
Buonanno (1985) also identified one such candidate star in For- 
nax 3. We note that these four stars (F3V01, F3V07, F3V88, and 
F3V94) appear centrally concentrated, with three having radial dis- 
tances r < 16". This suggests that at least some of these stars are 
cluster members, and worthy of further investigation. 

We can measure the intrinsic red and blue edges of the in- 
stability strip in each cluster from the colours of its RR Lyrae 
stars, and the colour excesses derived above. The four clusters seem 
particularly suitable for this type of measurement, with each pos- 
sessing a very well populated horizontal branch (excepting clus- 
ter 1). Examining each edge, we see that in general there is some 
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Figure 5. Colour-magnitude diagrams for each of the four clusters, showing their horizontal branch regions. The RR Lyrae stars are plotted. RRab stars are 
solid circles, and RRc stars solid triangles. Those stars with uncertain classification have hollow points. The candidate variable stars are plotted with crosses. 
Small dots represent non-variable objects. 



small mixing of variable and non-variable stars. In addition, in 
some cases (e.g., the red edge for Fornax 2) one variable star lies 
quite distinct from the remainder. To account for these two effects, 
we choose an "edge" to be the mean colour of the two most ex- 
treme variable stars. Using this definition, we have calculated the 
colours (V — I) re and (V — I)be of the red and blue edges for 
each cluster^. These quantities are listed in Table 2, along with 
their de-reddened counterparts. There is good agreement between 
the de-reddened measurements of the red edge. We obtain a mean 
(V — I)re = 0.59 ± 0.02. The blue edge is more poorly defined. 
The HB for Fornax 1 does not extend to non- variable stars, so we 
can discount the measurement for this cluster. Similarly, while the 
Fornax 2 HB has a significant population of BHB stars, there is a 
gap at the blue edge of the RR Lyrae strip. Hence, (V — I) be for 
this cluster is also ill-defined. For the well populated clusters For- 
nax 3 and 5, we obtain good agreement in the blue edge colour, 
with both having (V — I) be = 0.28. Applied to Fornax 2, this 
value has excellent consistency, with the reddest non-variable BHB 



1f We did not include the four "bright" variable stars in the calculation for 
Fornax 3. 



stars having intrinsic colours V — I ~ 0.27. We therefore adopt a 
mean (V - I) be = 0.28 ± 0.02. 

Finally, we can calculate the index (B ~ R)/{B + V + R) 
of Lee, Demarque & Zinn (1994) (hereafter referred to as the HB 
type). In this quantity, B is the number of BHB stars in a cluster, 
R is the number of RHB stars, and V the number of HB variable 
stars - in our notation these are Nbhb, Nrhb, and Nvhb re- 
spectively. In general, Nvhb Nrr, because we include N can d 
in the count, and have the possibility of excluding some of the stars 
counted in Nrr - e.g., the four bright variables in Fornax 3. If 
we initially ignore the possibility of field star contamination, we 
count Nrhb as the number of stars with 20.8 < V < 21.7 and 
0.59 < (V — I)o < 0.77. The red edge was chosen after careful 
consideration of the CMDs from the cluster centres, so as not to in- 
clude any RGB stars. The value is in excellent agreement with that 
adopted for M54 by Layden & Sarajedini (2000). The BHB stars 
are easier to count, since this is an isolated area in each cluster's 
CMD. We estimated the errors in Nrhb and Nbhb by count- 
ing the number of stars lying close to the edges of the defined re- 
gions. This gives some idea of the number of stars which might 
have been scattered in or out of these regions. We estimated the 
error in Nvhb by counting the number of non- variable stars ly- 
ing in the RR Lyrae strip, and adding this to N ca „d- This accounts 
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Figure 6. HB type versus metallicity diagram for the galactic globular clus- 
ters and the Fornax clusters from the present paper. The galactic clusters 
are split into three subsystems according to the criteria of Zinn (1993). 
The Fornax clusters are most similar to the galactic young halo popula- 
tion. The overplotted isochrones are from Rey et al. (2001) and constitute 
the latest versions of the Lee et al. (1994) synthetic HB models. The two 
lower isochrones are respectively 1 . 1 Gyr and 2.2 Gyr younger than the top 
isochrone. 



for stars which we may have mis-identified as variable, or variable 
stars which were missed by our detection procedure. All star counts 
are listed in Table 2. 

Using these quantities, we measured HB types of —0.30 ± 
0.04, 0.42±0.05, 0.40±0.05, and 0.52±0.04 for the four clusters. 
The errors were derived using the quadrature formula in Smith et 
al. (1996). Our measurements are in reasonable agreement with the 
values derived by Buonanno et al. (1998b) using the same WFPC2 
observations, but different numbers of variable stars (see Section 
3.3). However, there is the possibility that our counts as described 
above have been contaminated by field stars, especially for clusters 
2 and 3, which lie against significant backgrounds (see e.g., De- 
mers et al. (1994)). In particular, these two clusters seem to have 
quite well populated red clumps on the RGB near the RHB, and 
it is therefore likely our measured HB types are too red. We have 
no nearby field observations, so to try and account for this we re- 
calculated the HB index for each cluster using only stars measured 
on the PC, where the cluster centres are imaged. In this case we 
derived HB types of -0.38 ± 0.06, 0.50 ± 0.05, 0.44 ± 0.05, and 
0.52 ± 0.05, respectively. The calculation for cluster 1 involves 
small numbers of stars, and we prefer the original value. For the 
other three clusters, it is clear that Fornax 5 has not suffered any 
contamination (as expected - this cluster lies well separated from 
the dwarf galaxy), but Fornax 2 and 3 possess some small field star 
contamination. For these two clusters we prefer the HB types as 
calculated from the PC, rather than the original values. 

Additionally, we note that the CMD for cluster 1 presented 
by Buonanno et al. (1996) apparently shows several HB stars bluer 
than V — I ~ 0.3, whereas these stars are not present in our data. 



Hence our derived HB type for this cluster may be too red. It is 
possible that the WFPC2 observations did not include the region of 
Fornax 1 containing the BHB stars; however we find this explana- 
tion unlikely since the images cover ~ 75 per cent of the cluster, 
including its core. Even so, adding several BHB stars to the calcula- 
tion does not produce a large effect. Four stars would move the HB 
type to ~ —0.1, still in reasonable agreement with our preferred 
value. 

These measurements confirm the existence of a "second pa- 
rameter" effect among the Fornax globular clusters, as highlighted 
previously by several authors (Smith et al. 1996; Buonanno et al. 
1998b). For example, clusters 1 and 3 have very similar metallic- 
ities (see Table 2), but extremely disparate HB morphologies. In 
Fig. 6 we plot a HB type versus metallicity diagram to highlight 
this result. Lee et al. (1994) have suggested that age differences 
can mostly account for second parameter variations in HB type, in 
that younger clusters possess redder horizontal branches at given 
[Fe/H]. Zinn (1993) split the galactic globular clusters by metal- 
licity and HB type into three distinct groups. The most metal rich 
clusters are associated with the bulge and thick-disk of the galaxy, 
and the most metal poor with the halo. The halo clusters are subdi- 
vided into "young" and "old" groups on the basis of their HB type. 
From Fig. 6 it is clear that the Fornax clusters are most similar to 
the galactic young halo clusters. This is compatible with the asser- 
tion of Zinn (1993) that the young halo clusters have likely been 
accreted by the galaxy through merger events with satellite dwarf 
galaxies. 

Buonanno et al. (1998b) have measured the ages of Fornax 
clusters 1, 2, 3, and 5 to be identical to within ±1 Gyr, and equally 
coeval with the metal-poor galactic globular clusters M68 and M92. 
The isochrones overplotted on Fig. 6 show that based on HB type, 
the Fornax clusters should be more than 1.1 Gyr younger than the 
galactic old halo clusters, and moreover should have an internal 
age spread greater than ~ 1.1 Gyr (in that Fornax 1 should be 
younger). If the Buonanno et al. (1998b) ages are correct, then it is 
possible that age is not the sole second parameter operating in the 
Fornax clusters. Bearing in mind the extreme difference highlighted 
in Section 4.1 between the specific frequencies of RR Lyrae stars 
in the Fornax and galactic globular clusters, and the way that the 
HB index is defined, it seems likely that the relative number of RR 
Lyrae stars in a cluster has a significant bearing on the measured 
HB type. It is certainly interesting to note that the four galactic 
globular clusters with Srr > 60 are all members of the galactic 
young halo population. 



4.6 Distance and structure of the Fornax dwarf 

Using the mean V magnitudes we measured for the RR Lyrae stars, 
it is also possible to determine a distance modulus to each cluster. 
To do this, we need to know the absolute magnitude My(RR) of 
the RR Lyrae stars. This quantity is a function of metallicity, so 
that more metal poor stars are intrinsically brighter. We adopt the 
relation from the review of Chaboyer (1999): 



M V (RR) = 0.23([Fc/H] + 1.6) +0.56 



(5) 



which includes Hipparcos results. This is the same relation used by 
Layden & Sarajedini (Layden & Sarajedini 2000) to determine the 
distance to M54 and the Sagittarius dwarf galaxy. We note that the 
relation has an uncertainty of ±0.12 mag in the zero-point term, 
reflecting the discrepancies which exist between the different tech- 
niques to determine Mv(RR). 
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Table 4. Globular cluster and RR Lyrae population characteristics. 



Fornax 1 



Fornax 2 



Fornax 3 



Fornax 5 



Nrr 

N cand 



15 
1 



43 
3 



99 
6 



40 

3 



M v 



(Lower limits) 



-5.33 
118.2 



-7.05 
70.5 



-7.80 
90.2 



-7.05 
68.4 



N„ 



N c /Nrr 

(P ab ) 

Ptr 
(Pc) 

«y-i) m ) 



E(V - /) 
E(B — V) 

A v 

(V - I) RE 

(V-I)be 

Nbhb 
Nrhb 

NvHB 
HB index 



[Fe/H] 
My (RR) 

<V(RR)> 



(m - M) 
Distance 



(All RRab) 

(RRab with good P) 

(RRab with AV > 0.9) 



(All RRab) 
(RRab with good P) 
(RRab with AV > 0.9) 
(All RRc) 

(All RRab) 

(RRab with good P) 

(RRab with AV > 0.9) 



Measured 

De-reddened 

Measured 

De-reddened 

(= B) (All chips) 

(= R) (All chips) 

(= V) (All chips) 

(i+V+fl) (All chips) 
(PC only) 



(All RR) 

(All RRab) 

(RRab with good P) 

(RRab with AV > 0.9) 

(All RRc) 



(kpc) 



10 
5 
4 
5 
0.33 

0.546 
0.611 
0.544 
0.431 

0.68 ±0.01 
0.70 ±0.01 
0.69 ±0.01 
0.10 ±0.01 
0.07 ±0.01 
0.23 ±0.02 

0.68 ±0.02 
0.58 ±0.03 
0.45 ±0.01 
0.35 ±0.02 

0±0 

7± 1 

16 ±2 

-0.30 ±0.04 
-0.38 ±0.06 

-2.05 ±0.10 
0.46 ±0.02 

21.27 ±0.01 
21.29 ±0.01 
21.29 ±0.01 
21.29 ±0.01 
21.26 ±0.01 

20.58 ±0.05 
130.6 ±3.0 



30 
15 
10 
13 
0.30 

0.494 
0.574 
0.540 
0.373 

0.65 ±0.01 
0.66 ±0.01 
0.66 ±0.01 
0.07 ±0.01 
0.05 ±0.01 
0.16 ±0.02 

0.66 ±0.02 
0.59 ±0.03 
0.39 ±0.02 
0.32 ±0.03 

68 ±4 

14 ±3 

46 ±9 

0.42 ±0.05 
0.50 ±0.06 

-1.83 ±0.07 
0.51 ±0.02 

21.34 ±0.01 
21.32 ±0.01 
21.31 ±0.01 
21.29 ±0.01 
21.36 ±0.01 

20.67 ±0.05 
136.1 ±3.1 



61 
32 
14 
38 
0.38 

0.532 
0.613 
0.582 
0.404 

0.63 ±0.01 
0.63 ±0.01 
0.63 ±0.01 
0.05 ±0.01 
0.04 ±0.01 
0.12 ±0.02 

0.66 ±0.02 
0.61 ±0.03 
0.33 ±0.02 
0.28 ±0.03 

114 ± 15 
20 ±4 

101 ± 15 

0.40 ±0.05 
0.44 ±0.05 

-2.04 ±0.07 
0.46 ±0.02 

21.24 ±0.01 

21.25 ±0.01 
21.21 ±0.01 

21.20 ±0.01 

21.21 ±0.01 

20.66 ±0.05 
135.5 ±3.1 



13 

8 
5 

27 
0.68 

0.532 
0.577 
0.548 
0.374 

0.62 ± 0.01 
0.60 ±0.01 
0.60 ±0.01 
0.04 ±0.01 
0.03 ± 0.01 
0.10 ±0.02 

0.62 ± 0.02 
0.58 ± 0.03 
0.32 ± 0.01 
0.28 ± 0.02 

83 ±6 

12 ±3 

43 ±5 

0.52 ±0.04 
0.52 ± 0.05 

-1.90 ±0.06 
0.49 ± 0.02 

21.33 ±0.01 
21.35 ±0.01 
21.32 ±0.01 
21.30 ±0.01 
21.30 ±0.01 

20.74 ±0.05 
140.6 ±3.2 



To obtain metallicities for each cluster, we used the recent 
measurements and literature collation of Strader et al. (2003). From 
their Table 5, we calculated the mean measured metallicity for each 
cluster. This combines the results of several different measurement 
techniques, ranging from integrated spectroscopy to the slope of 
the RGB. For the four clusters, we obtain [Fe/H] = -2.05 ± 0.10, 
-1.83 ± 0.07, -2.04 ± 0.07, and -1.90 ± 0.06. This results 
in RR Lyrae absolute magnitudes of My (RR) = 0.46 ± 0.02, 
0.51 ± 0.02, 0.46 ± 0.02, and 0.49 ± 0.02, respectively. 

The next step is to obtain a mean V RR Lyrae magnitude, 
(V(RR)) for each cluster. We determined these from the (V) mea- 
surements in Table 2. Stars with (V) deviant by more that 2<r were 
not included in the final average. The mean V magnitudes so ob- 



tained are listed in Table 4. Using the extinction values calculated 
above, together with the RR Lyrae absolute magnitudes, we then 
determined a distance modulus for each cluster. These results are 
also listed in Table 4. The mean distance modulus is 20.66 ± 0.03 
mag, which is in good agreement with previously determined val- 
ues. In particular, Buonanno et al. (1999) obtain an estimate of 
20.68 ± 0.20 to the centre of the Fornax dwarf (i.e., cluster 4) after 
a thorough discussion of several techniques and literature measure- 
ments. We note that the uncertainties quoted for our results are due 
to random errors only, and that there are also systematic errors to- 
talling approximately ±0.15 mag. 

In addition we briefly investigated the likelihood that our 
(V(RR)) values are systematically too faint, based on the discus- 
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sion in Section 3.4. We recalculated (V(RR)) using only RRab 
stars, then only RRab stars with good periods, and then only RRab 
stars with A V > 0.9. These values are also listed in Table 4. On av- 
erage, the (V(RR)) calculated using all the RRab stars are ~ 0.03 
mag fainter than those calculated using the RRab stars with good 
periods, and ~ 0.04 mag fainter than those calculated using the 
RRab stars with AV > 0.9. This demonstrates exactly what we ar- 
gued in Section 3.4 - that we expect a systematic error in (V(RR)) 
due to our short observation baselines. If we include the RRc stars 
in our calculations, the offsets are moderated somewhat. The two 
worst affected clusters are Fornax 3 and 5, which may have distance 
moduli 0.03 mag closer than those we have quoted. 

Our distance moduli correspond to distances of 130.6 ± 3.0 
kpc, 136.1 ±3.1 kpc, 135.5 ±3.1 kpc, and 140.6 ±3.2 kpcforthe 
four clusters. Again, the quoted uncertainties reflect random errors 
only. These measurements are consistent with, but do not require, a 
significant line of sight depth to the Fornax system, with a value of 
~ 10 kpc. The outlying clusters are Fornax 1 and 5, which are also 
the clusters with the greatest angular separation from the centre of 
the Fornax dwarf (see Mackey & Gilmore (2003b)). We note that 
if Fornax 5 is indeed 0.03 mag closer than we have listed above, 
the line of sight depth shrinks slightly to ~ 8 kpc. Demers et al. 
(1994) obtained a tidal radius for Fornax of r t — 77 ± 10 arcmin 
along the minor axis, and an ellipticity of 0.27 at this radius. Hence 
the tidal radius along the major axis must be 105 ± 14 arcmin. 
These radii correspond to linear distances of 7.0 ± 0.9 kpc and 
9.6 ± 1.3 kpc at our mean distance modulus of 20.66. These are of 
a similar magnitude to the ~ 4.5 kpc radius obtained from our line 
of sight depth, suggesting that Fornax is close to spherical. This 
lack of extended line of sight depth argues against models of dSph 
galaxies which invoke very large line of sight depth and/or tides to 
explain the observed high stellar velocity dispersions. Other recent 
studies have reached a similar conclusion (e.g., Klessen, Grebel & 
Harbeck (2003), who studied the thickness of the Draco dSph HB). 
It seems most likely that dSph galaxies with high stellar velocity 
dispersions constitute dark matter dominated systems. 



RR Lyrae stars. The specific frequencies of the four clusters are 
exceptionally large in comparison to those for the galactic globu- 
lar clusters. All of the four Fornax clusters have Srr > 80, while 
only four of the ~ 150 galactic clusters have such large values. It 
is likely that Fornax 1 has the largest specific frequency ever mea- 
sured in a globular cluster. The Fornax clusters are also unusual 
in that their RR Lyrae populations possess mean properties which 
lie intermediate between the two Oosterhoff groups defined by the 
galactic globular clusters. In this respect the RR Lyrae stars in the 
Fornax clusters are similar to those in the old LMC clusters, but 
most resemble the field populations of several dwarf galaxies. For- 
nax cluster 5 stands out in having an extremely high fraction of 
RRc stars (N c /Nrr ~ 0.6), and possessing four strong RRe (sec- 
ond overtone pulsator) candidates. 

Using the mean RR Lyrae colours at minimum light, we have 
measured colour excesses and line of sight extinction values to- 
wards the four clusters. With a significant catalogue of HB variable 
stars, we have also been able to provide new estimates of the HB 
morphology in the these clusters. We have confirmed the presence 
of the second parameter effect, in particular between clusters 1 and 
3, which have identical metallicities and ages, but extremely differ- 
ent HB morphology. The HB morphologies of the Fornax clusters 
are in general redder than many galactic clusters of similar metal- 
licity - in this respect they most resemble the "young" galactic halo 
population. We also find good agreement between measurements of 
the intrinsic red and blue edges to the instability strip at the HB. We 
determine {V-I) RE = 0.59±0.02 and (V-I)be = 0.28±0.02. 

Finally, we have used the mean RR Lyrae brightnesses to de- 
termine the distance to the Fornax dwarf galaxy. We measure a 
mean distance modulus of (m — M)o = 20.66 ± 0.03 mag, which 
is in excellent agreement with previous measurements. Our calcu- 
lations are consistent with a line of sight depth of ~ 8 — 10 kpc for 
the Fornax dwarf galaxy. This value is in good accordance with the 
dimensions of this galaxy as measured in the plane of the sky, and 
argues against the tidal-remnant models of dSph galaxies which are 
invoked to explain their observed internal stellar dynamics. 



5 SUMMARY AND CONCLUSIONS 

We have identified and measured 197 RR Lyrae stars in four of the 
globular clusters in the Fornax dwarf galaxy, using archival WFPC2 
observations. We also located 13 variable HB objects, which are 
likely RR Lyrae stars. A typical star in our sample has 14 paired 
F555W and F814W observations spread over ~ 0.36 days. Such 
a short baseline has compromised our ability to measure accurate 
periods for ~ 50 per cent of the RRab type stars we identified - 
for these objects we instead calculated lower limits for the oscil- 
lation periods and amplitudes. Nonetheless, we have successfully 
determined periods for all the RRc and 50 per cent of the RRab 
stars in the ensemble. We also present measurements of mean mag- 
nitudes and colours for all the objects in our sample, and discuss 
the introduction of systematic errors into these calculations due 
to our short observational baseline. This is the first time that RR 
Lyrae stars have been identified in these clusters, and quantitative 
measurements of their variability presented. Naturally, additional 
observations at more than one epoch will be extremely useful in 
refining our classifications and measurements, and in particular in 
allowing good periods to be determined for all stars. 

Even given the difficulties associated with measuring periods 
for some of the sample, we have been able to obtain a large amount 
of useful information about the Fornax globular clusters from the 
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